Nonadiabatic Channels in the Superconducting Pairing of Fullerides 
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We show the intrinsic inconsistency of the conventional phonon mediated theory of superconduc- 
tivity in relation to the observed properties of RbaCeo- The recent, highly accurate measurement 
of the carbon isotope coefficient ac = 0.2i, together with the high value of Tc (30 K) and the very 
small Fermi energy Ef (0.25 eV), unavoidably implies the opening of nonadiabatic channels in the 
superconducting pairing. We estimate these effects and show that they are actually the key elements 
for the high value of Tc in these materials compared to the very low values of graphite intercalation 
compounds. 

PACS number(s): 74.20.-z, 74.7Q.Wz, 7f.38.+i 



One of the most striking evidences of the phonon role in 
high-temperature superconductivity of alkali-doped Cgo 
compounds is the observation of nonzero isotope effect 
on the value of the critical temperature Tc [Q. However, 
the large spread of the reported values of the carbon iso- 
tope coefficient ac = —d\aTc/d\aM where M is the 
isotopic mass, has prevented to settle a definitive and 
self-consistent picture. This long standing uncertainty 
has been solved only recently for the compound RbaCeo 
(Tc ~ 30 K). Resistive measurements on 99% enriched 
-'^^C single- crystals have permitted in fact to establish 
ac =0.21 with high precision |^. 

The knowledge of the accurate value of ac is an im- 
portant element to establish not only the important 
role of phonons, but also to test the self-consistency of 
the Migdal-Eliashbcrg (ME) theory of superconductiv- 
ity in alkali-doped fuUcrenes. The measured val- 
ues ac = 0.21 and Tc = 30 K can be used to extract 
the microscopic quantities involved in the superconduct- 
ing pairing. For example, in Ref. |3|, Tc = 30 K and 
ac = 0.21 are interpreted within the conventional ME 
theory by A = 0.9, win = 1360 K and ^x* = 0.22, where A 
is the electron-phonon coupling constant, win is the log- 
arithmic phonon frequency and fi* is the Coulomb pseu- 
dopotential According to this standard analysis, the 
high value of Tc in alkali-doped Cgo compounds is merely 
due to a strong electron-phonon coupling to the high- 
est intramolecular phonon modes. These results should 
be compared with the graphite intercalation compounds 
(GIC) where Tc ~ 0.2 K |7| is explained by a moder- 
ate coupling (A ~ 0.3) to similar high energy phonon 
modes. Current theories claim that the big difference 
between the electron-phonon coupling in fullerides com- 
pared with graphite intercalation compounds arises from 
the finite curvature of the Ceo molecule . In this per- 
spective, therefore, RbsCgo is just an ordinary strong- 
coupling superconductor described by the conventional 
adiabatic ME framework. 

In this Letter instead we demonstrate the intrinsic in- 



consistency of the standard ME theory in RbsCgo. This 
conventional description is in fact invalidated by the ex- 
tremely low value of the Fermi energy Ep — 0.25 eV 
~ 2900 K characteristic of the fuUerene compounds [Q. 
We show that the whole range of A-Wph values which 
fit Tc = 30 K and ac = 0.21 through the solution of 
the standard ME equations implies a Migdal parameter 
Xujph/Ep larger than 0.4, instead of being zero as as- 
sumed by the ME theory . This situation inevitably 
leads to the breakdown of Migdal's theorem and to the 
opening of nonadiabatic channels in the superconduct- 
ing pairing By solving the nonadiabatic equations 
||l^,0, we estimate these effects and show that they are 
actually the key elements for the high values of Tc in the 
fuUerene compounds. 

We now discuss why the experimental data Tc = 30 K 
and ac = 0.21 ^ are inconsistent with respect to the ME 
theory. To simulate the coupling of the electrons to the 
different intramolecular phonon modes, we consider an 
electron-phonon spectral function a'^F{uj) modeled as a 
rectangle of width Aluq and centered at ojq. For Awq = 0, 
a'^F{uj) reduces to a single Einstein (5-peak while for 
Awo > it becomes broadened. The electron-phonon 
coupling constant is determined by the usual relation A = 
2 J duja'^F{uj)/uj. The Coulomb pseudopotential is taken 
to have the standard form /i* — /i/[l -t- /i ln([jjc/wniax)]j 
where Wmax = ^^o + AciJo/2 is the maximum phonon fre- 
quency, fj, is the screened Coulomb parameter and uJc 
is the high-frequency cut-off. According to whether it 
is the entire set of 7r-band s E2| or rather only the nar- 
row conducting tiu band ||l^f \vhich contributes to the 
dynamical screening, we shall consider lUc = 5wmax or 
UJc = E-p = 2900K, respectively. For different values of 
Awo/wq we then solve numerically the ME equations to 
find the values of A, luq, and fj, (or fi*) which reproduce 
the experimental data Tc = 30 K and ac = 0.21. 

In Fig. |] we show the calculated fi and /i* (top panel) 
and ujq (lower panel) as function of A for Aluq/luq — 
(solid lines) and Awq/wq = 1 (dashed lines) with ujc — 
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FIG. 1. Plot of the Coulomb parameters n and (up- 
per panel) and of the central phonon frequency tJo (lower 
panel) obtained by the standard ME equations constrained 
to have Tc — 30 K and ac = 0.21. Solid lines: AtJo/<^o = 0, 
LOc = 5(Jn,ax; dashed lines: Acjq/^^o = 1, tJc = Su^maxl 
dot-dashed lines: Aljq/ujo = 0, u>c ~ Ep. 



Swmax- The main point of Fig. ^ is that the calcu- 
lated depends strongly on the elcctron-phonon con- 
stant A. For large values of A, Tc — 30 K and ac = 0.21 
are reproduced only for quite small phonon frequencies 
while decreasing A quickly enhances luq. The Cgo phonon 
spectrum however is limited by a maximum frequency of 
~ 2300 K @Jl|l and this settles a lower limit for the al- 
lowed values of A. In fact, from the lower panel of Fig. 
|l|, it is clear that for A smaller than 1 the corresponding 
^max = ^0 + Acjq/S rapidly exceeds 2300 K signalling 
that the solution of the ME equations falls well outside 
the range of applicability for the fuUerene compounds. 

A further interesting point is that the effect of the spec- 
tral broadening (Awo > 0) is of secondary importance 
to the overall behavior, suggesting that the results are 
only weakly affected by the detailed structure of a'^F{uj). 
In fact, for A = 0.9 we find )i* ~ 0.25, wi„ ~ 1873 
K (Awq = 0) and cjin — 1555 K (Awq = "^o), where 
Win — exp[(2/A) / duj In uja'^F{Lo)/uj] is the logarithmic 
phonon frequency which, for the rectangular model here 
considered, is given by cjin = [ujq — (AtJo/2)^]^^^- These 
values are consistent with those reported in Ref. ||] 
(A ~ 0.9, fi* ~ 0.22, Win — 1360 K) where a phonon spec- 
trum obtained by ab initio calculations has been used to 
fit Tc = 30 K and ac = 0.21. 

In Fig. 1^ we show also the results for Acjq — and uic — 
Ep (dot-dashed lines). Due to the reduced dynamical 
screening of the Coulomb repulsion, the increase of uiq is 
much steeper than the previous cases until at A 1.25 we 
find ojq ~ Ep and /i* ~ /i so that it is no longer possible 
to have ac smaller than the BCS value 0.5. 

We now address the consistency of the data of Fig. 
|l| with the standard ME theory [Q|5j. The assumption 
at the basis of the ME framework is Migdal's theorem 
which states that, as long as the phonons have a much 
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FIG. 2. Values of LOp^i/Ep (upper panel) and of 
P — \upi-i/Ep (lower panel) extracted from the data of Fig. 
1 with Ef = 0.25 eV= 2900 K. Solid lines: Einstein phonon 
spectrum; dashed lines: broad spectrum with Awo/ij-'o = 1, 
dot-dashed lines: Einstein phonon spectrum Awo/wo = 1 
with LUc = -Ep. The Migdal-Eliashberg theory holds true only 
when ujph/Ep <g 1 and P <^1. 



slower dynamics than that of the electrons, the nonadi- 
abatic electron-phonon interference effects (vertex pro- 
cesses) can be neglected This condition is well sat- 
isfied in conventional superconductors since their Fermi 
energy is of order Ep ^ 10 eV ^ 10^ K while the highest 
phonon frequencies are usually less than ^ 50 meV ^ 60 
K 1^. However, the alkali-doped fuUerene compounds 
are molecular solids characterized by very narrow con- 
duction electron bands of width of only W ~ 0.5 eV 
||l|. In RbaCeo (as in the other AaCeo compounds) the 
conduction band is half-filled by electrons and the corre- 
sponding Fermi energy is Ep ~ 0.25 eV — 2900 K, while 
the maximum phonon frequency is ~ 2300 K. In princi- 
ple, therefore, there is no reason to expect Migdal's the- 
orem to be applicable in fuUerene compounds, unless the 
main interaction is with the lowest Cgo phonon modes 
('^ 400 K) via however a rather weak coupling. 

We can test whether the data of Fig. |l are consistent 
with Migdal's theorem by evaluating the order of magni- 
tude (P) of the first nonadiabatic electron-phonon vertex 
correction. By following the same reasonings of Migdal 
|0,|o|, this is given by: 



P = 2 / dhj !^ 

Jo Ep 



A 



(1) 



where for the rectangular model we considered the aver- 
age phonon frequency cjph reduces to 



duj a^F[Lo) 



Awo 



ln( "" + ^"°/g ) 



(2) 



When P ^ the nonadiabatic interferences are negligi- 
ble and the ME theory holds true; on the contrary, size- 
able values of P signal the breakdown of the standard 
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theory. In Fig. || we show the values of LOp^/Ep (upper 
panel) and of P (lower panel) extracted from the ujq val- 
ues reported in Fig. || and by setting Ep — 2900 K. As 
expected from the overall trend of ujq vs A plotted in Fig. 
|l|, the adiabatic ratio uj-ph/ Ep is large and close to unity 
for A < 1, while it rapidly decreases to ujph/Ep ~ 0.1 
for large values of A. One could therefore argue that 
adiabaticity is guaranteed for very large electron-phonon 
couplings and that in this regime the ME framework is 
valid. This is however incorrect because according to 
the vertex correction is proportional to A so that, as 
shown in the lower panel of Fig. ^ P is never negligible. 
Note that the claimed value A ~ 0.9 1^ corresponds to 
ujpii/Ep > 0.7 and a minimum P > 0.6, i.e., the vertex 
correction is comparable to unity. 

From Figs. |l| and ^ we conclude therefore that the 
conventional phonon-mediated superconductivity is not a 
complete and self-consistent picture of RbsCeo since the 
values of A and wq needed to fit Tc — 30 K and ac — 0.21 
strongly violate Migdal's theorem. This conclusion holds 
true even when electron-phonon spectra more structured 
than the rectangular one are used to fit the data of Ref. 

. By adding a 5-peak centered at w = to a rectangu- 
lar spectrum of width covering the whole intramolecular 
modes, we have in fact simulated additional contributions 
from low-frequency (wi ~ 50 K) Cgo-Ceo phonon modes 
[ p5[ , and from an enhanced coupling to soft intramolec- 
ular modes ||T^ {coi ~ 400 — 600 K) possibly related to 
dynamical Jahn- Teller effects 0. We find that = 30 
K and ac = 0.21 imply P > 0.4 when tJi = 50 K and 
P > 0.45 when cui — 400 K (further details will be pre- 
sented elsewhere). 

The above results point out that, if superconductivity 
in RbaCgo is mediated by phonons, a consistent descrip- 
tion of its superconducting properties should be sought 
beyond the ME theory. More precisely, the low value of 
Ep indicates that the adiabatic hypothesis and Migdal's 
theorem should be abandoned from the start and that the 
theory should be formulated by allowing ujpi^/Ep to have 
values sensibly larger than zero. This naturally leads to 
nonadiabatic interference effects in the electron-phonon 
scattering which can significantly modify both the nor- 
mal and superconducting properties with respect to the 
ME phenomenology 

Indeed, characteristic effects of the nonadiabatic ver- 
tex corrections are predicted to be observable in several 
quantities, like Tc and its the isotope coefficient P,|l0|, 
the reduction rate of Tc itself upon disorder ||ll| , the ef- 
fective electronic mass m* |lj], the Pauli susceptibility 
pgf . A peculiar feature of the nonadiabatic processes is 
to produce, under some conditions, constructive electron- 
phonon interference in the particle-particle channel lead- 
ing to an enhancement of Tc [p|JlC|]. Hence Tc — 30 K, 
for a given phonon spectrum, can be achieved by much 
smaller values of A than needed in conventional ME the- 
ory. Favourable conditions to this trend are expected in 
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FIG. 3. Plot of /I, /I* and uio calculated by the nonadiabatic 
theory as solutions of Tc = 30 K and ac ~ 0.21. Note that 
the range of values for A is now much smaller and realistic 
with respect to those of Fig. |l|. 

materials with strong electronic correlation, as fuUerenes: 
strong local repulsion suppresses short-range interactions 
(large q's in Fourier space) and favours forward small-q 
scattering | |20| , |2l[ | for which electron-phonon vertex pro- 
cesses become attractive. 

Now we re-analyze the experimental constraints of 
RbaCgo, Tc = 30 K and ac — 0.21, in the context of the 
nonadiabatic theory of superconductivity. We show that 
the inconsistencies of the results derived by ME theory 
are naturally solved when the same experimental data 
are coherently analyzed in the nonadiabatic regime. 

Explicit analytical and diagrammatic equations of the 
nonadiabatic theory of superconductivity have been out- 
lined in previous works H-|l^ and, for sake of shortness, 
they will be here omitted. A set of generalized Eliash- 
berg equations in nonadiabatic regime is constructed by 
following a perturbative approach based on the Migdal 
parameter P. The consistency of such a perturbative 
scheme is discussed below. We simplify the phonon spec- 
trum by assuming a dispersionless Einstein phonon with 
energy wq- We note however that our results, as also 
shown in Figs. and |[ are only weakly affected by 
the specific shape of a'^F{uj). Electronic correlation is 
taken into account by a cut-off Qc on the electron-phonon 
exchanged momenta, which selects forward scattering, 
where stronger the correlation smaller qc. 

In Fig. ^ we show the phonon frequency luq, the stat- 
ically and dinamically screened Coulomb repulsion, re- 
spectively and /i*, vs. A obtained in the nonadiabatic 
theory to reproduce Tc = 30 K and ac = 0.21. The 
parameter qc has been chosen qc = 0.2/eF (fcp the Fermi 
vector), an appropriate value for a strongly correlated 
system, and the dinamical screening of fj, has been con- 
sidered to be provided by the tiu electrons H. 

From the comparison of Fig. || with Fig. ]l|7a first re- 
markable difference lies in the distinct ranges of electron- 
phonon couplings needed to reproduce the RbaCeo data. 
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FIG. 4. Adiabatic ratio ujp]^/ Ef and vertex correction mag 
nitude P obtained by the nonadiabatic solutions of Fig. ^ 



In fact, the conventional ME theory predicts A Jc^ 1 (Fig. 
|l|), while the nonadiabatic analysis yields A 1 (Fig. ||). 
But the most striking difference is that, if we now take 
the parameters obtained by the generalized theory and 
use them in the standard ME theory, these would give 
a very low value of Tc, less than 1 K or even zero. This 
result is now perfectly compatible with the GIC super- 
conductors, for which Migdal's theorem holds true, and 
it clarifies that the high Tc values of the fuUerides are 
essentially due to constructive nonadiabatic interference 
effects rather than to a very large value of A. In our per- 
spective, therefore, the origin of the enhancement from 
Tc ~ 0.2 K in GIC to ~ 20 - 30 K in fuUerene com- 
pounds stems mainly from the opening the nonadiabatic 
channels in the electron-phonon interaction, rather than 
from a ^ 300% enhancement of A. 

We address now the consistency of the perturbative 
scheme with the nonadiabatic solutions for RbaCgo- In 
Fig. ^ we show the adiabatic ratio cUph/Ep {ujph = ojq) 
and the Migdal's parameter P extracted from the data re- 
ported in Fig. ^. The large value of ajpii/Ep shown in the 
upper panel points out again the breakdown of Migdal's 
theorem and consequently the need of the inclusion of 
the nonadiabatic vertex corrections. The magnitude of 
the vertex corrections P ^ 0.2, certainly not negligible, 
is however small enough to support a perturbative ap- 
proach in P §,0. Note moreover that, according to the 
comparison with exact results for the single-electron Hol- 
stein model p^ , for weak couplings the system is away 
from polaron formation and that the perturbative scheme 
is well defined. 

In conclusion, we have investigated the validity of 
Migdal-Eliashberg theory of superconductivity in RbaCgo 
by analyzing the constraints imposed by recent experi- 
mental data, namely the critical temperature Tc — 30 
K and the isotope effect ac = 0.21. We have found 
that the values of A and Wph needed to reproduce the 
experimental data, together with the very low value of 



the Fermi energy, strongly violate Migdal's theorem and 
are therefore inconsistent with the ME framework. This 
situation unavoidably leads to the opening of nonadia- 
batic channels in the electron-phonon pairing which we 
argue to play the primary role for the high values of Tc in 
fuUerene compounds. Finally, we stress the importance 
of peculiar nonadiabatic effects in both superconducting 
| pT[ and normal state properties ^8 19 1 of fuUerene com- 
pounds. Experiments in this direction are therefore of 
great interest. 
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